Brahmagupta's interpolation formula